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Abstract. I present detailed models for the formation of disk galaxies, and inves- 
tigate which observables are best suited as virial mass estimators. Contrary to naive 
expectations, the luminosities and circular velocities of disk galaxies are extremely poor 
indicators of total virial mass. Instead, I show that the product of disk scale length and 
rotation velocity squared yields a much more robust estimate. Finally, I show how this 
estimator may be used to put limits on the efficiencies of cooling and feedback during 
the process of galaxy formation. 



1 Introduction 

Currently, the main uncertainties in our picture of galaxy formation are related 
to the intricate processes of cooling, star formation, and feedback. The cooling 
and feedback efficiencies are ultimately responsible for setting the galaxy mass 
fractions /gai — -Mgai/Afvir- Here Mgai is the total baryonic mass of the galaxy 
(stars plus gas, excluding the hot gas in the halo) and A/vir is the total virial 
mass. 

Here I present new models for the formation of disk galaxies, which I use to 
investigate how well observables extracted directly from these models can be used 
to recover /gai(-Afvir). Even though the assumptions underlying the model are not 
necessarily correct, and the phenomenological descriptions of star formation and 
feedback are certainly oversimplified, this provides useful insights regarding the 
ability of actual observations to constrain the poorly understood astrophysical 
processes of galaxy formation. 

2 Short Description of Models 

The main assumptions that characterize the framework of the models are the fol- 
lowing: (i) dark matter halos around disk galaxies grow by the smooth accretion 
of mass, (ii) in the absence of cooling the baryons have the same distribution of 
mass and angular momentum as the dark matter, and (iii) the baryons conserve 
their specific angular momentum when they cool. I follow Firmani & Avila- Reese 
(2000) and make the additional assumptions that (iv) the spin parameter of a 
given galaxy is constant with time, (v) each mass shell that virializes is in solid 
body rotation, and (vi) the rotation axes of all shells are aligned. Although nei- 
ther of these assumptions is necessarily accurate, it was shown in van den Bosch 
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(2001) that they result in halo angular momentum profiles in excellent agreement 

with the high rcsohition iV-body simulations of Bullock at al. (2001). 

The main outline of the models is as follows. I set up a radial grid between 
r = and the present day virial radius of the model galaxy and follow the 
formation and evolution of the disk galaxy using a few hundred time steps. Six 
mass components are considered: dark matter, hot gas, disk mass (both in stars 
and in cold gas), bulge mass, and mass ejected by outflows from the disk. The 
dark matter, hot gas, and bulge mass are assumed to be distributed in spherical 
shells, whereas the disk stars and cold gas are assumed to be in infinitesimally 
thin annuli. Each time step the changes in these various mass components in 
each radial bin are computed using the following prescriptions. 

• The rate at which the total virial mass grows with time is given by the 
Universal mass accretion history (van den Bosch 2002a). 

• At each redshift, the dark matter is assumed to follow an NFW density 
distribution (Navarro, Frcnk & White 1997). 

• The gas that enters the virial radius of a halo is added to the disk a time 
tc = max[fff,tcooi] later, where tg and tcooi correspond to the free-fall time 
and cooling time. The cooling time depends on the metallicity of the hot 
gas, which is taken to be a free parameter Zhot- 

• The radius at which the gas settles is governed by its specific angular mo- 
mentum distribution, which follows from the assumptions (iv)-(vi) listed 
above. 

• When the disk becomes unstable, part of the disk material is converted into 

a bulge component (cf. van den Bosch 1998). 

• In the disk, only the cold gas with a surface density above the critical density 
given by Toomre's stability criterion is considered eligible for star formation. 
This gas is transformed into stars with a rate given by a simple Schmidt law. 

• Part of the cold gas is expelled from the system by supernovae feedback. 
This is modeled as in van den Bosch (2000), and regulated by a free feedback 
efficiency parameter efb- 

More details regarding these models can be found in van den Bosch (2001), 
where it is shown that these models yield disk galaxies in good agreement with 
observations. For instance, for the majority of the model galaxies the disk re- 
veals an exponential surface brightness profile. Note that contrary to previous 
disk formation models (e.g.. Mo, Mao & White 1998), this is not an a priori 
assumption of the model. 

3 Galeixy Mass Fractions 

In order to investigate how /gai(-Mvir) relates to the cooling and feedback effi- 
ciencies 1 discuss three models that only differ in the metallicity of the hot gas, 
Zhot, and the feedback efficiency, efb: the 'Standard Model' with Zhot = 0.3 Zq 
(a typical value for the hot gas in clusters) and £fb = (i.e., no feedback), the 
'Zero Metalhcity Model' with Zhot = 0.0 and £fb = 0, and the 'Feedback Model' 
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with Zhot = 0.3^0 and £fb = 0.02 (i.e., two percent of the SN energy is converted 

to kinetic energy). All other model parameters are kept fixed at their fiducial 
values (see van den Bosch 2001). For each model a sample of 400 model galaxies 
is constructed. Present day virial masses are drawn from the Press-Schechter 
mass function with 10^"h~^ Mq < Mvir(O) < lO^^/i"^ Mq, corresponding to 
31 kms~^ < Kir < 312 kms~^, roughly the range expected for galaxies. Spin 
parameters, which parameterize the specific angular momentum, are drawn from 
a typical log-normal distribution. 

Figure 1 plots the present day galaxy mass fractions /gai = Mgai/Myir as 
function of Mvir (crosses). In the models without feedback /gai is virtually iden- 
tical to the universal baryon fraction /bar for low mass systems. For more massive 
systems, cooling becomes ineflScient, causing /gai to strongly decrease with in- 
creasing virial mass. This is more pronounced in the model with Zhot = 0, for 
which cooling is least efficient. In the Feedback Model (right panel), /gai ^ /bar 
for the low mass systems, but with a large amount of scatter. This is a reflection 
of the scatter in halo spin parameters: systems with less angular momentum 
produce disks with higher surface densities, therefore have higher star formation 
rates, which induce a more efl[icient feedback. At the high mass end /gai is fairly 
similar to the Standard Model without feedback. This owes to the fact that the 
mass ejection efficiency scales inversely with the square of the escape velocity, 
making feedback less efficient in more massive systems. 
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Fig. 1. A comparison of the true galaxy mass fraction /gai as function of the true virial 
mass (crosses) with the same values estimated from the observables extracted from 
the models (dots). Results are plotted for all three models discussed in the text. The 
horizontal dotted line corresponds to the universal baryonic mass fraction /bar . In the 
model with feedback (right panel) the dots occupy the same parameter space as the 
crosses, indicating that tfie observables allow one to recover /gai(Afvir), at least in a 
statistical sense. In the two models without feedback there are significant errors in the 
recovered /gai(A/vir). Yet, the two /gai(Afvir) are sufficiently different to discriminate 
between the two models; in particular, the estimated galaxy mass fractions nicely avoid 
the upper right regions of parameter space which contain information on the cooling 
efficiencies. Furthermore, the dots occupy different areas of parameter space in models 
with and without feedback, such that there is hope that the "observed" /gai(Mvir) may 
be used to constrain the efficiency of feedback. 
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Clearly /gai(Afvir) depends strongly on both the cooling and feedback effi- 
ciencies. Therefore, if one could obtain a measure of /gai(-Mvir) observationally 
it would allow us to constrain the poorly understood physics of cooling and 
feedback. This requires one to be able to infer both the baryonic galaxy mass 
Mgai as well as the total virial mass Afvir from observations of the luminous (and 
gaseous) components. Using the models outlined above I now investigate which 
observables are best suited as the appropriate mass indicators. 

For the galaxy mass one can write Mgai = TbLb + Mgas. Here Lb is the 
total i3-band luminosity, Tg is the corresponding stellar mass-to-light ratio, and 
Mgas is the galaxy's (cold) gas mass. is easily obtained observationally, and 
I assume that Mgas is observationally accessible through HI measurements. The 
stellar mass-to-light ratio, however, is not directly accessible to an observer, but 
has been shown to correlate strongly with color. I therefore extract the B — K 
color from the models from which I estimate Tb using the relation of Bell & de 
Jong (2001). 

The galaxy virial mass is more difficult to obtain. In Figure 2 I plot three 
'observables' as function of Afvir for galaxies in each of the three models. Note 
that both Vmax, defined as the maximum rotation velocity inside the radial extent 
probed by the cold gas, and Lk are poor indicators of virial mass. First of all 
the slope and zero-points of the Knax(Mvir) and LxiM^ir) relations depend on 
the input parameters of the model. This means that an observer trying to infer 
Mvir from either Vinax or Lk needs to make assumptions about the efficiencies 
of cooling and feedback. However, it is exactly these efficiencies that we seek to 
constrain. Secondly, the scatter of both relations can be so large that even if the 
normalization of the relation were known, one could still not infer Mvir to better 
than an order of magnitude. In particular, the scatter in Vinax(-Afvir) can be very 
large. This owes entirely to the scatter in halo spin parameter, which sets the 
concentration of the baryonic mass component after cooling, therewith strongly 
influencing T/max- 

As the lower panels in Figure 2 show, a much more reliable virial mass in- 
dicator is RdV^^^/G. Here Rd is the disk scale length in the /-band, obtained 
from fitting an exponential to the /-band surface brightness distribution of the 
disk. Upon fitting all galaxies of all three models simultaneously I obtain 

M.„ = 2.54xlO.»M,(-|.)(^J|^)\ (1) 

(rms scatter between 20 and 50 percent, depending on the amount of feedback). 
The fraction of model galaxies for which equation (|l|) yields an estimate of 
the true virial mass to better than a factor two is larger than 97 percent! It 
is remarkable that the zero-point for models with feedback is the same as for 
models without feedback. When matter is ejected it reduces Vmax but at the same 
time increases the disk scale length such that /JdVj^ax stays roughly constant. 
This is due to the star formation threshold criterion included in the models. If 
real galaxies follow a similar threshold criterion, equation (1) provides a fairly 
accurate estimate of the total virial mass of disk galaxies. 
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We now have all the tools in place to see whether we can recover /gai(Afvir) 
from the "observables" . The dots in Figure 1 correspond to the estimates of 
/gai(Mvir) obtained using the method outlined above. In both models without 
feedback there are significant errors in the recovered /gai(Afvir), which is domi- 
nated by errors in the estimate of Mvir- Yet, the recovered /gai(Afvir) of the two 
models are sufficiently different to distinguish between them. In particular, in 
both models the dots avoid the regions in the upper right corner which contains 
information about the efficiency of cooling. In the feedback model the recovered 
values occupy roughly the same area of the /gai — M^ir plane as the intrinsic 
values. Although the one-to-one correspondence for individual model galaxies 
may be poor, statistically the method to recover /gai(Afvir) explored here works 
reasonably well. 

4 Conclusions 

We have shown that the product of disk scale length and maximum rotation 
velocity squared can be used as a fairly accurate estimator of total virial mass. 
This in turn can be used to obtain estimates of the galaxy mass fractions as 
function of virial mass, which contains important information about the effi- 
ciencies of cooling and feedback. Another approach, that has been taken in the 
past, is to use published luminosity functions and luminosity- velocity relations 
to construct halo velocity functions (i.e., Newman & Davis 2000; Gonzales et 
al. 2000; Kochanek 2001). The main goal of these studies is similar to the work 
presented here, namely to circumvent the problems with poorly understood as- 
trophysical processes when linking the observed properties of galaxies to those 
of their dark matter halos. Our results imply that great care is to be taken in 
linking an observable velocity such as Vmax to the circular velocity of a dark 
matter halo. Based on our results, we suggest that the construction of a halo 
mass function using A/vir oc i?ciKnax may proof more reliable. More details of 
the results presented here can be found in van den Bosch (2002b) 



References 



9. 

10 
11 
12 
13 



1. E. F. Bell, R. S. de Jong: ApJ, 550, 212 (2001) 

2. J. S. Bullock et al.: ApJ, 555, 240 (2001) 

3. C. Firmani, V. Avila-Reese: MNRAS, 315, 457 (2000) 

4. A. H. Gonzales et al.: Ap J 528. 145 f2000) 

5. C. S. Kochanek: preprint, |astro-ph/0108160| (2001) 

6. H. J. Mo, S. Mao, S. D. M. White: MNRAS, 295, 319 (1998) 

7. J. F. Navarro, C. S. Frenk, S. D. W. White: ApJ, 490, 493 (1997) 

8. J. A. Newman, M. Davis: ApJ, 534, Lll (2000) 
F. C. van den Bosch: ApJ, 507, 601 (1998) 
F. C. van den Bosch: ApJ, 530, 177 (2000) 
F. C. van den Bosch: MNRAS, 327, 133 4 (2001) 



F. C. van den Bosch: MNRAS, in press; ^stro-ph/010515^ (2 002a) 
F. C. van den Bosch: MNRAS, submitted~|astro-ph/0lT2566| (2002b) 




Fig. 2. The relation between various virial mass estimators and the actual virial mass 
for all three models. In the upper panels log(Knax) is plotted versus log(Mvir)- The 

1/3 

dotted line corresponds to Mvir oc Knax, which is a reasonable description of the av- 
erage relation. However, the scatter is large, and the zero-point depends on the actual 
model, which makes Knax unsuitable as virial mass indicator. The same goes for the 
K-hand luminosity, which is plotted in the middle row of panels. Here the dotted line 
corresponds to Lk tx Afvir, which only yields a reasonable description for the brighter 
galaxies in models with .Zhot ~ Zq /3. For fainter galaxies, however, the scatter is large. 
Furthermore the slope of the Lit (K/ir) relation depends strongly on the feedback effi- 
ciency. Thus total luminosity is also a poor indicator of total virial mass. The lower 
panels plot the virial mass estimator RdV^^^/G as function of Mvir. Here the dotted 
lines correspond to equation (|l|), which gives a reasonable fit to the model galaxies, 
independent of the cooling and/or feedback efficiencies. In addition, the scatter is rel- 
atively small, making the product of disk scale length and maximum rotation velocity 
squared a fairly accurate estimator of total virial mass. 



